Context
This project aims at exploring how various metrics/proxies can give insights on the vertical distribution of the phytoplankton biomass in the water column using the data measured during the Green Edge oceanographic cruise conducted in the Baffin Bay in 2016. Four main indicators were considered to get insights about the phytoplankton biomass in the water column:
- Fluorescence
- Pigments
- Particle beam attenuation coefficient (CP)
- Particle backscattering coefficient (bbp)
Most of the data is presented as a function of open water days (OWD) that were calculated in Randelhoff et al. (2019). The first graph shows the OWD of each station.
A note on the CTD data
After exploring the CTD data, I found out that there were some unwanted peaks in both the florescence and transmittance data. Therefore, I have used a median moving window to smooth the data (both transmittance and fluorescence).
Data processing for 3D plots
During the meeting we had on 2021-01-22, it was questioned why I was averaging data by owd and depth before doing the interpolation. This is simply because it is possible that there are many observations for a specific pair of owd and depth. For example, it can be seen here that there are 9 values (i.e. 9 stations) that have an owd = 2 and depth = 1.980.
| owd | depth_m | N |
|---|---|---|
| 2 | 1.980 | 9 |
| 2 | 2.178 | 9 |
| 2 | 2.376 | 9 |
| 2 | 2.574 | 9 |
| 2 | 2.772 | 9 |
| 2 | 2.970 | 9 |
Here are 9 stations in the CTD data with owd = 2 and depth = 1.980.
Hence, before performing the interpolation, data has to be averaged by OWD and depth.
Fluorescence
This section shows the vertical distribution of phytoplankton based on fluorescence.
Fluorescence from the CTD
Using boxplots can be one interesting alternative way to present the data. For example, I have divided the CTD fluorescence data presented in the above graphic into ice-free/ice-covered and above/below the 0.1 isolume.
Fluorescence from the MVP
Note that because the MVP is measuring continually, there are no stations associated with each measurement. This is why that the MVP data is not presented as a function of OWD.
Pigments
These graphs were made using the pigments data from the rosette. Pigments were summed into two groups:
- Photoprotection pigments
- Accessory pigments
Particle beam attenuation coefficient (CP)
The increase of cp during the lighted portion of the day is usually explained by the accumulation of intracellular carbon concentration associated with photosynthetic processes (Kheireddine and Antoine 2014).
CP from the CTD
The CP visualization based on the CTD can be seen above. There is a clear relationship between fluorescence and CP as seen below.
The next graph shows the same data but divided differently. The best relation can be seen for the observations below the 0.1 isolume and located in open water. I think the bottom left panel is the most interesting where we can see the L shape.
Finally, the same data by transect.
The next graphic shows CP vertical profiles for some of the deepest stations below 500 meters. Any ideas why there are almost always a bump in the data in the very last portion of the profiles?
CP from the MVP
The next graphs show the longitudinal variability of CP at the deepest locations within each transect. The left column shows how the MVP was moving in the water column whereas the right column shows how CP varied along the transect. For example, we can see that the MVP followed the 300 meters line at transect 100.
Particle backscattering coefficient (bbp)
bbp from the hydroscat
These graphs show bbp at six different wavelengths measured by the hydroscat.
Particulate organic carbon (POC)
bbp vs chlorophyll-a fluorescence (Hydroscat)
The next graphs explore the relationships between bbp and chlorophyll-a fluorescence from the Hydroscat device.
Same data divided differently.
Comparing bbp at two wavelengths (Hydroscat)
Scatterplot between bbp measured at 532 nm and 700 nm.
Same data presented differently.
Phytoplankton absorption
Maybe we could use absorption information to get some additional insights.
Vertical profiles
This graph shows the vertical profiles of phytoplankton absorption at 440 nm. It seems that open water stations have subsurface maximum more apparent than the profiles of the ice-covered stations.
Spectral profiles
By looking at the spectral profiles, we can also notice that there are indeed differences in the water column.
Finally, I have averaged the spectral profiles in four categories. I think these can be seen as end-members spectra.
Comparing blue vs red regions
With all the data.
By groups.
blue/red ratio of specific absorption
We can see that there is a an increase of the ratio after owd = 0. I am not sure to understand why the ratio is also higher at depth when owd < 0. Is there a shift when the ice is melting?
POC vs CP
The next graphs show the relationships between POC and CP.
bbp vs chla (Hydroscat)
Note that I have used chla fluorescence from the Hydroscat which is given as raw counts. If we decide to go further with this data, I guess we could convert it in actual biomass/stock quantities.
bbp vs cp
The ratio bbp(532)/cp(660) could be considered as a proxy of parti- cle size and composition, increasing when small or inorganic particles become relatively more abundant than large or organic particles (Xing et al. 2014).
So, based on the next graphs, larger particles are more abundant at depth than at the surface. However, I do not really see trends in the data.
xxx
This index was calculated as folow:
\[ \frac{\log(\frac{\text{bbp}(532)}{\text{bbp}(700)})}{0.274} \]
